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Abstract 

We study a two-dimensional charged particle interacting with a mag- 
netic field, in general non- homogeneous, perpendicular to the plane, a 
confining potential, and a point interaction. If the latter moves adia- 
batically along a loop the state corresponding to an isolated eigenvalue 
acquires a Berry phase. We derive an expression for it and evaluate it 
in several examples such as a homogeneous field, a magnetic whisker, 
a particle confined at a ring or in quantum dots, a parabolic and a 
zero-range one. We also discuss the behavior of the lowest Landau 
level in this setting obtaining an explicit example of the Wilczek-Zee 
phase for an infinitely degenerated eigenvalue. 



1 Introduction 



A nontrivial Berry phase | Bei ] can be demonstrated in different situations. 
There is a growing interest recently to this effect in mesoscopic systems - see, 



e.g., IILSQ , |MHK|| and references therein. These papers investigate theoreti 



cally and experimentally how the phase is manifested in quantum dynamics 
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of a particle with spin interacting with a time-dependent magnetic field. In 
the present paper we are going to discuss a simple model in which the Berry 
phase emerges even if the spin-orbital coupling is neglected. 

The model describes a charged particle confined to a potential well and 
placed into a magnetic field of constant direction, which is independent of 
time and may be homogeneous. The phase will appear when the well is mov- 
ing adiabatically. The similar situation appears in the Born-Oppenheimer 
approach for the study of molecules (see, e.g., and references therein) 
and impurities in semiconductors ||Za2|| . For the sake of simplicity we suppose 
that the well represents a zero-range interaction, i.e., it is given by a point 
interaction in the plane. This makes it possible to derive explicit formulae for 
the Berry potential. The idea of employing point interactions to this purpose 
is not new: some solvable models exhibiting a nontrivial Berry phase have 
been constructed earlier. For instance, the geometric phase resulting from 
a cyclical motion of the boundary condition for the Dirac and Schrodinger 
equations on an interval [0,^] was computed in ||BFG| , |GK|| . On the other 
hand, the authors of ||CS|| investigated the Berry phase which arises when a 
pointlike scatterer is adiabatically moved in a rectangular billiard in such a 
way that the energy levels encircle a "diabolic point" . 

In our cases the results are simpler and rather illustrative. In particular, 
we shall show that moving the zero-radius potential well along a closed curve 
in the plane, the eigenfunctions of a particle trapped by the well and exposed 
to a homogeneous field perpendicular to the plane acquire a phase which co- 
incides with the number of magnetic field quanta through the area restricted 
by the curve. This picture changes if an additional confining potential is 
added, say, in the form of an annular potential "ditch" . In the limiting case 
of an infinitely thin ring the motion of the point interaction induces a ge- 
ometric phase which differs from the above one on a quantity proportional 
to the persistent current in the annulus. Recall that persistent currents in a 
ring with a point perturbation were investigated - see, e.g., [ |(JGK| - but the 
relation to the Berry phase was not noticed . 

Let us describe briefly the contents of the paper. In the next section 
we shall recall briefly how the zero-range interaction in a magnetic system 
is constructed and how its spectrum is determined by means of the Krein's 
formula. For simplicity we suppose always that the magnetic fleld as well 
as the possible conflning scalar potential are rotationally symmetric. The 
central part of the paper is Section 3 where we derive a general expression 
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for the Berry potential corresponding to a point interaction moving along a 
smooth curve - cf. Eq. (p.l8|) . 

This result is in the next section illustrated on the number of examples. 
We show that the Berry phase for the perturbation moving along a closed 
loop C in a homogeneous field without a scalar potential is proportional to 
the number of flux quanta through C. In distinction to that, the phase cor- 
responding to a magnetic whisker contains an extra term proportional to the 
persistent current in the loop C. For comparison we analyze an electron con- 
fined to a circular ring and find the same Berry phase expression containing 
the persistent-current term, in this case independently of the field profile. Fi- 
nally, we discuss a harmonic quantum dot in a homogeneous field. We show 
that if the point interaction is strong the effect of the confining potential is 
small and the Berry phase is again given by the number of the flux quanta 
through C, up to an error term. We compare this with the situation where 
the quantum dot itself is zero-range. 

The behaviour of degenerate eigenvalues under adiabatic change of pa- 
rameters is more complicated and less understood. In the final section ad- 
dress this question in the present setting and discuss what happens in this 
situation with the lowest Landau level. We compute the generalized Berry 
potential which determines the corresponding Wilczek-Zee phase, and find 
the latter for adiabatic evolution around a small loop. It appears to be 
nontrivial for the angular momentum m = 1 while the states with higher 
momenta are not affected. Moreover, the phase which arises here differs in 
sign from the one corresponding to the isolated energy level; we explain this 
effect ClS Sb sort of topological charge conservation. 

2 Magnetic systems with a point 
perturbation 

As indicated above we shall consider a charged particle of charge e and mass 
m^, (which may be thought of as the effective mass of an electron in a crystal) 
living in the plane with Cartesian coordinates x, y and exposed to a magnetic 
field perpendicular to the plane, B = B{x,y)ez. We also assume that the 
particle may be confined to a part of the plane by a non-negative potential 
W. The main simplifying assumption we shall make concerns the rotational 
symmetry: we suppose that there is a system of polar coordinates r, (p such 
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that the magnetic field and the confining potential depend on the radial co- 
ordinate only, B = B{r) and W = W{r). In this case one can choose a 
gauge in such a way that the radial component of the vector potential van- 
ishes, Ar{r,ip) = 0, and A^{r,ip) = A^{r) depends on r only. In particular, 
WA = 0. 

It is convenient to single out the uniform component of the magnetic field, 
B = Bq + Bi with Bo being a fixed vector. Of course, such a decomposition 
is arbitrary, but we will have mostly in mind situations when B has a finite 
limit as r — oo; then the non- uniqueness is removed by the requirement 
Bi ^ 0. We shall also employ the corresponding decomposition of the vector 
potential, A = Aq + Ai. In view of the assumed symmetry it is natural to 
use the circular gauge, Aq^v) = ^B^r e^. As for the nonconstant part, we are 
particularly interested in the example of an infinitely thin Aharonov-Bohm 
solenoid, or a magnetic fiux line with y4i^(r) = where $ is the magnetic 
fiux through the solenoid. It is convenient to use a dimensionless parameter 
T], T] = $/$0) where 

27[hc 

\e\ 

is the magnetic fiux quantum; so rj is the number of quanta carried by the 
solenoid. The corresponding magnetic field is concentrated at the origin of 
the coordinates, Bi = <l>5(r)e^. 

In the following considerations, however, we assume only that Ai^p is a 
smooth function of the r variable on the halfiine (0,oo). Under the stated 
assumptions, the particle Hamiltonian has the following form: 



e 
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H = -ihV - -A] +W 

2m* V c / 



2 



2m^, 
2m* 



1 d d 1 
■r— — h 



+ (b + -A (r)] (t— 

r dr dr ' r"^ dcp"^ \ 2m*c \ r ^'^ J \ dip 



(^^Sy + BorAi^ir) + Ai^{ 



ry]+W{r). (2.1) 



The potential W is non-negative by assumption, and therefore if is a well 
defined self-adjoint operator, which can be understood, e.g., as the Friedrichs 



extension of the operator defined on C^(]R^ \ {0}) by the rhs of Eq. (2.1). 

Now we shall introduce a point perturbation of the above Hamiltonian 
located at a point s G with the polar coordinates (p, 9). (Further we will 
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assume that sj^ if Ai has a singularity at the point r = 0). The perturbed 
operator H^^g is conventionally obtained as a self-adjoint extension of the 
symmetric operator S which is a restriction of H to the domain 

V:={^eV{H): ^{s) = 0}; (2.2) 

since the deficiency indices of 5* are (1, 1) the extensions are characterized by a 
single parameter a. Under rather general assumptions about the Hamiltonian 
GMQI the Green function Ga,s{r, r'; E) of Ha,g is given by the Krein formula 



G^,,{f, f E) = G{f, f E) - [Q{E; s) + a]-'G{f, s; E)G{s, f E) , (2.3) 

where Q{E] s) is the so-called Krein Q-function or renormalized Green func- 
tion at the diagonal point [s, s), 

Q(E;s) := lim G(r,s;E) - — ^ln|r- , (2.4) 

and a is the mentioned parameter. The latter is related to the scattering 
length A of the point interaction by the formula 

a = !Ih_^\n\-\ (2.5) 



^2 



Less formally, the point perturbation at a point s may be defined via the 
Fermi pseudopotential of the form 

/x(5(f - s) (1 - In |f - s |(f - s)V?) , (2.6) 

where the coupling constant /i is related to the parameter a by /i = a~^. 

Under rather weak regularity requirements on the potentials A{r) and 
W{r) the Green function is of the form 

rn 

G{f, f'; E) = — ^ In |f - r"|-i + G'o(f, f'; E) , (2.7) 

where Go is continuous in the variables r, r ' and analytic with respect to E in 
the resolvent set, C \ <y{H), of the free operator. It is the case, for example, 
if Ai and W are smooth functions (see, e.g., [p^j . Chapter III, Theorem 
5.1.). If Ai has a singularity at the origin, then every point r, r 7^ 0, has a 
neighborhood such that Eq. ( |2.7|) is true for f' in this neighborhood. 
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s, s 



E) 



(2.8) 



Since the singular term in Eq. ( |2.7| ) is energy independent, dG/dE needs 
no renormalization and we have 

dQ{E-s) _ dG 
dE dE 

Due to the well-known Weyl theorem, the essential spectra of H and H^^g 
coincide. As for the discrete spectrum, it may happen that H and Ha^g have 
a common eigenvalue. Let E be an isolated eigenvalue of H such that there 
exists a corresponding eigenfunction if) satisfying ip{s) = (in particular, 
this can be always achieved if is a degenerate eigenvalue). Then E belongs 
to the spectrum of Ha^g as an isolated eigenvalue; moreover, the multiplicity 
m' of E in the spectrum of Ha^g obeys the inequality m' > m — 1, where m 
is the multiplicity of E in the spectrum of H. This assertion may be proven 
following the arguments from ||CdV|| where a special case of our claim has 
been considered. In addition, the spectrum of Ha,g contains all solutions of 
the equation 

Q{E;s) + a = (2.9) 

(the true levels of the zero-range well). Every solution of this equation lies 
in a gap of the unperturbed spectrum and is a simple isolated eigenvalue of 
Ha^g. The corresponding eigenfunction ip has the form 

-1/2 



Gif,s;E) 



(2.10) 



Recall that in the real part of the resolvent set the derivative is positive, 
{dQ/dE){E) > for E G M \ a{H) - cf. Therefore, equation 

has at most one solution in every gap of the spectrum cr{H). Generally 
speaking, the equation (|2.9| ) may have no solutions - see, e.g., ||AGM|| . It is 
straightforward to see that if Eq is an isolated eigenvalue of H and ip{s) ^ 
holds for at least one eigenfunction corresponding to Eq, then Eq is a pole of 
the function Q(-; s). Hence if <7{H) is purely discrete solutions of Eq. ( |2.10| ) 
exist in infinitely many spectral gaps. 

A simple but important particular case of the considered problem, Ai = 
and W = 0, concerns a free motion in a uniform magnetic field. In this 
situation the Green function acquires the following explicit form. 



G{f,f';E) 



2^ 



E 



exp 



-iri^Qr A r' 



[r — r 



/\2 



E 

huJr 



(2.11) 
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where ujr 



UJr_ 



is the cyclotronic frequency, 



le^ol 



eBo 
27ihc 



is the flux density of the uniform component of the magnetic field, ao. 



h 



-1/2 



(2.12) 



is the magnetic length, and \1' is the Tricomi confluent hypergeometric func- 
tion - c/. [ pMM|| . The Q-function now obviously does not depend on s and 
equals |GHS| : 



Q{E) 



2^ 



E 

hjjJr 



+ 27 - ln2 - 21nao 



(2.13) 



where 'il){x) = (lnr(a;))' and 7 = —ipil) is the Euler constant. Up to a 
scaling and a shift in the argument the behaviour of Q{E) is given by that 
of the digamma function ip; this shows that in a uniform magnetic field the 
zero-range potential with any fixed a G M induces existence of an energy 
level on the halfline [—00,60) as well as in each interval (e^,5^+i), where 
Ee := + I) are the Landau levels. 



3 The Berry phase 

Let us return to the condition (|2.9| ). In what follows we will keep a fixed 
(and drop it mostly from the notation) and move the point s* along a smooth 
path C : s = s(t),t G [0,1], in the plane (or in the punctured plane 
\ {0} if Ai has a singularity at the point 0) in such a way that Eq. ( |2.9|) 
has a solution Eo{s,a) lying in a gap of the unperturbed Hamiltonian H. 
Denote 

-1/2 

G{f,s;Eo{s,a)). (3.1) 



dQ 
dE 
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the corresponding normalized eigenfunction of the perturbed operator Ho,^s 
(see (|2.10|) . If the path C is a closed loop, s (0) = s (1), the initial and final 

state, 'ips{Q) and ips{i), respectively, differ by a phase factor, 

#(1) = ^s(o) exp {-'- ^ E^{s{t)) dt + 27(C)^ , (3.2) 

where the Berry phase 7(C) depends only on the path C; in accordance with 
Ref. [[Berj it equals 

7(C) = / V{s) ds , (3.3) 
Jc 

where 

V{s) :=z(#| V.-lV's), (3.4) 

is the so-called Berry vector potential. Recall that from the differential- 
geometric point of view Im {ipg \ | ips) is a connection 1-form in a principal 
fiber bundle over (or M? \ {0}) associated with the eigenfunction fibration 
ipg^ s [0]; in other words, this quantity is a gauge potential with the gauge 
group U(l). We shall express V{s) in the polar coordinates, 

V{s) = Vp{p, 9)e, + Ve{p, 9)ee (3.5) 

with 

Vp = I Vp I ijg) , Ve='- I Ve \ • (3.6) 

To proceed further we need more information about the structure of the 
Green function G{f, r'\E). First of all, we decompose the state space L^(]R^) 
into partial waves, i.e., we represent it as L^(M+,rdr) ® L^(§^,d(y9) and 
perform the Fourier transform on the second component, L^(S^, d^p) ^^(Z) 
with 



9'^{gm}m&, gm = ^= g{(p) 6 '"""^dip . (3.7) 

V 27r Jo 

Then L^(]R^) decomposes into an orthogonal sum of subspaces each of which 
is isomorphic to the radial component, 

oo 

L2(m2) ^ L2(M+,rdr). (3.8) 



m=— oo 
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The unperturbed operator H commutes with rotations around the origin, 
and therefore it decomposes correspondingly into the orthogonal sum 



H 



(3.9) 



where the partial-wave parts Hm are self-adjoint operators in L^(M+,rdr) 
obtained as the Friedrichs extensions of the operators (|2.1|) with the domain 
C^(M+, r dr) and —id/ dip replaced by m. It is obvious that each Hm is a real 
operator, i.e., that it commutes with the operator of complex conjugation in 
L^(M+, r dr) . It follows that its Green function Gmi'f', r'; E) is real valued for 
a real E. 

The full Green functions can be expressed through its partial-wave com- 
ponents as 



G{f,f';E) = — e^"^('^-'^'^G^(r,r';E). 



(3.10) 



m=— oo 



It follows that {ipg I Vp I ips) is a real number. On the other hand, differenti- 
ating the identity {ipg \ i'g) = 1 we see that the real part of {ipg \ Vp | tjjg) (as 
well as (ipglVdlipg)) vanishes. Consequently, the radial component of the 
Berry potential Vp{s) = 0. To find the angular one, let us differentiate the 
identity 

{H-E)G{f,s;E) = 6{f- s) (3.11) 

with respect to the constant component Bq of the magnetic field keeping E 
and s fixed; this yields 



0. 



Notice that dG/dBo is a smooth function in view of ( p.7| ). Hence 

dG 



(3.12) 



G 



(H-E) 



dG{f, s; E) 



dBo 



dBo 
dG 



dBo 



s, s; E) 



dQ 
dBo 



{E, s \ 



and therefore 



G 



dH 
dBo 



G 



dQ 
dBo 



(3.13) 
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Dividing both terms of this expression by OQ/dE and putting E = Eq{s), 
we arrive at the relation 



OH 

dB'o 



dQ fog 



dBo \dE 



(3.14) 



Since Eq{s) solves the equation (PT^), the last term at the Ihs can be expressed 

as 



dQ fdQ 



dBn \dE 



so 



dH 
dB'o 



E=Eo 



OEq 
dBo' 



dEo 
dBo 



Now we shall employ the formula ( |2.1| ) which yields 
OH . eh d e^Bo 2 



dBo 2m^:C dip ~'~ Am^:C 



rAi^(r) 



(3.15) 



(3.16) 



^ G{f, s; E), and furthermore, 



It follows from (^) that ^G{f, s; E) = 
that Q{E,s) is independent of the angular variable, Q{E,s) = Q{E, p). The 
last claim means that Eo{s) also does not depend on 9. As a result we have 
that Veips = — V^Vs- Finally, we express the angular momentum operator 
from ( |3.16| ) as 



. d dH ^ ^ , , 27r ^ , , 

dip eh oBo $0 

which allows us to cast the sought angular component into the form 
1 



(3.17) 



P 



2m^c dEois) ^ , , , 2 1 , \ / \ 2vr , , . ^ / s , , > 
+ 7r^o(^s I r I ips) + (sgn e) — {i)s \ rAi^{r) \ ips) 

<I>o 



eh dBo 
dEo{s 



271 



Tih^ d^o 



+ vr^o(V's \r \ipg) + (sgn e) — {ipg \ rAi^{r) \ ipg) 



$0 



We stress that in view of (|3.10D Ve is independent of 9. 



(3.18) 



4 Examples 

Let us now illustrate the Berry phase behaviour on several examples. 
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4.1 A homogeneous field 

Suppose that the magnetic field is uniform, i.e., Ai = 0. Since the Green and 
Krein functions are exphcitly known in this case, it is convenient to evaluate 



the Berry potential directly from the relation ( |3.4| ). It follows from (|2.11| ) 
that ipg is of the form 

ijg{f) = exp [-7rz^o(r A s)] /(|f - s\) , (4.1) 

and therefore 

^ei^si^) = — niC,Q rp {cos ip cos 9 + sin ip sin 9) exp [—niC,Q{r A s)] f{\r — s\) 
+ exp [-TTiUrA s)] Ve/(|r - s\) . (4.2) 

Inspecting the explicit form of the function / we see that it is real-valued 
and normalized, = 1. It follows that 

(^,1 exp l-mUrA s)] Ve\f) = {f\ V, | /) = ^Ve||/f = , 
so the sought quantity is given by the first term only. 



{ijs\Ve\iJs) = -7rt^o / f-s\f{\f-s\)\'df 

= -7iz^o f {s^ + r-s)\fi\r\)\'df (4.3) 



= — vri^ofp / 1/(1^1)1 df+p / r(cos(y9cos^^ + siny9sin6')|/(|r 1)1 dr). 

The first integral obviously equals one and the second zero, hence 

Ve{p) = -{^s\Ve\^s)=7r^oP (4.4) 
P 

and the Berry phase is given by 

7(C)=27reo5, (4.5) 
where S is the area encircled by the loop C. We can write it also as 

7(C) = 27rsgne^, (4.6) 
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where is the full magnetic flux through the loop. Comparing ( [4 .41 ) which 
corresponds to Ai = with the general expression (|3.18|) derived in the 
previous section, we get in the limit p ^ the relation 



^c(^c?k1^o)- (4.7) 



Let us finish the example with a remark concerning an extension of the above 
result to three-dimensional systems. Suppose that the field is parallel to the 
z-axis and B{p,(p,z) = B{p) holds in the cylindrical coordinates. Then 
we have Vp = 0, Vq = vr^oP? and = 0, where (p, 9, () are the cylindrical 
coordinates of the point s. Consequently, the Berry phase along a closed loop 
C is again 7(C) = 2n^c/^o, up to a sign, where is now the magnetic 
flux through the projection of C to a plane perpendicular to the field. 



4.2 A magnetic whisker 



The opposite extreme corresponds to the situation where the homogeneous 
component is absent and the field is concentrated into a flux line (sometimes 



called Aharonov-Bohm solenoid), i.e 
yields 



Bo = and = ^ e^. Then ( CT) 



V 



Ve{p) = -7^ + - sgne . 



(4.8) 



Suppose, in particular that the point perturbation moves along a circle C of 
radius R centered at the origin of coordinates. In that case the Berry phase 
equals 

2m* dEo 



7(C) = - 



+ 2tt (sgn e) t] . 



(4.9) 



Bo=0 



The total flux of the field B through the circle C is {nR^^o + r])^o. 

a 



Keeping the flux $ fixed, we have ^ — --e)2^_ 8 



nR^^QjSr-- Hence 



7(C) 



7. 7^ \- 27r sgn e 

^2 d^c $0 



(4.10) 



Bo=0 



Recall that for a particle confined to the loop C the derivative dEo/d^c 
equals — ^/o where Iq is the corresponding persistent current. To understand 
better the meaning of the Eq. ([4.10|) we consider the following example. 
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4.3 Electron in a ring 



Up to now the confining potential of (2J.) was trivial. Tlie previous example 
inspires us to analyze another extreme situation in which is a very deep 
and narrow well. To get a solvable model we employ the usual idealization 
and suppose that the particle is confined to an infinitely thin circular ring 
C pierced by the magnetic field. In that case the Hamiltonian H becomes 
one-dimensional. Having in mind an electron, e < 0, we can write H as 



H 



. d 

dip 



(4.11) 



where R is the ring radius and $ = ?7$o is the total flux of the field B through 
the circle; the field profile is irrelevant here. The Green and Krein function 
are of the form 



-,i'm{ip—ip') 



(m + rjy — 



E 



(4.12) 



and 



m^R 



m=— oo 



2m^R^ 
f9 



E 



-1 



sin^V2m7E 



hy/m^E cos ^^y/2m^E — cos 2TTr] 
Consider now a point perturbation of the operator H, 

He = H + a-^6{ip - 6) . 
As above the Green function for Hq is given by the Krein formula 
Ge(<^, <^'; E) = G(<^, <^'; E) - [Q{E) + a]-^G{ip, 9; E)G{9, <^'; E) . 
A solution to the spectral condition 

Q{E) + a = 



(4.13) 

(4.14) 
(4.15) 
(4.16) 



exists in each interval {Ei, Ee^i), where {Ei}i>Q is the sequence of "free" 
eigenvalues 
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arranged in the ascending order. In addition, for a < the Eq. ( 4.1(j| ) has a 
solution also on the halfiine {—oo,Eq). 

Consider a fixed solution Eo{6) of ( [4.16D . It is clearly independent of 9 
and represents a nondegenerate eigenvalue of Hq with the eigenfunction 



dQ 
dE 



.-1/2 



G{^,d-E, 



(4.18) 



Let us evaluate the Berry phase when the perturbation site 9 travels once 
around the ring. The Berry potential is given by 



V{9) = i{M'^e\i^e) 



We express ipe in the form 



[m + rf) 



2 2m, 771 



(4.19) 



(4.20) 



with CO := [§{Eo)Y'^' 
V{9) 



Then 



On the other hand, 
dQ 



dE 



(Eo) 



7l¥ 



so 



- E 



m=— oo 



E 



((m + r]) 



(4.21) 



(m + 77) 



2 2Tn^R 



2 1-2 



(4.22) 



E 



m=— 00 



. m=— 00 



((m + r/)2 - ^^E) 
Differentiating now Q with respect to r] we get 



H -2 



E 



-1 



2m, 



9?7 



E 



m + ?7 



m=— 00 



(4.23) 



(4.24) 
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which yields the identity 

oo 

m 



E 



nfi^ dQ 



2m, 



2m^R dr] 



f)r} ' J 



V 



2 • 



In combination with (|4.22|) and ( [4.23|) this formula gives 



V 



(4.25) 



(4.26) 



h? dr] 

and the corresponding Berry phase accumulated while 6 moves once anti- 
clockwise around C is 



7(C) 



(4.27) 



h? dr] 

Taking into account that the total flux through the ring is = ^7^0 we see 
that the obtained expression is fully analogous to the formula ( [4.10| ) valid in 
the whisker case. 



4.4 A parabolic quantum dot 



As the next example of this section we shall discuss a quantum dot in a 
uniform magnetic fleld Bq. To get a solvable model, we suppose that the 
conflning potential which determines the dot is parabolic, W{r) 
The frequency ujq is related to the effective radius R of the dot by 



2^2 



c 



-m^:UjQR^ 



(4.28) 



where C, is the chemical potential of the system ||BL|| . The spectrum of H is 
discrete with the eigenvalues (usually called the Fock-Darwin levels) 



\m\ + 1 



+ n ) + hujcU , m G Z, n G N - 



(4.29) 



where u := y^cu|~+u^. We can employ the known propagator kernel | |K(J| | of 
the operator H, 



K{f, r'; t) 



m^uj 



rm^uj 



inih sin ^ \ ih sin ^ l 

—IT ■ r cos lir A r sm 

2 2 



/ 2 I i2\ 

(r + r ) cos 



ijjt 



(4.30) 
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To find an integral representation of tlie Green functions of H, one has to 
perform the Wick rotation in ( |4.30| ), i.e., to pass to the imaginary time t — >■ 
—it. This yields the heat kernel of e~*^; applying the Laplace transformation 
to it we get 



rriM 



27rh? Jq ^ [ Ahsinhiot 

-2r ■ r' cosh ujct + 2if A r ' sinh ujct 



+ r cosha;t 



sinh ujt 



(4.31) 



We shall also need the Krein function. It is obtained by the following trick: 
we observe that replacing uoc at the rhs of ( 4.31 ) by uj we get the Green 
function of the Landau Hamiltonian with the cyclotronic frequency uj. We 
add and subtract this function at the rhs , then we subtract the singularity, 
In 1?^ — r '|~^, and pass to the limit f*, r ' — > s. In accordance with (|2.13|) 
we obtain 

Q{E-s) 



2^ JO 
dt 

X 



^2tE/h 



sinh cut 2Txffi 



2h sinh cot 
1 E \ 
2~h^)^ 



{cosh out — coshtUc^) — 1 



27-ln2-21na 



(4.32) 



where a :- 



h 



We shall not analyze the last expression generally and restrict ourselves 
to showing that if the point-interaction is strong enough in the sense that 
Eq <^ —huj, the confinement potential has an insignificant effect on the Berry 
potential only. To this aim we denote 2E/h = —e and split the integral / 
in rhs of ( [4.31| ) into a sum I = Ii{e) + /2(^) of integrals corresponding to 
the intervals (0,e~^/^) and {e'^^"^, oo). It is easy to see that the first integral 
obeys the inequality hie) > ci{f,f')e~^ with a constant Ci depending on 
r and r ' only. On the other hand, using an integration by parts we find 
that hi^) < C2(r, r")^~^e^^. Neglecting for large \E\ the second integral, 
we have 



G(r,r ;rj ^ tttt? ^^P 



27rh^ 



-I—. — r A r 



2h 



X 



-1/2 



e^^^/'^exp 



[r — r 



l\2 



dt 

T 



(4.33) 
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Since the integral depends on 



/|2 



only, we can repeat the considerations 



of Section 4.1 obtaining thus 



Va 



T:i,p + 0{\E\o')- 



(4.34) 



4.5 A zero-range quantum dot 



The results of the previous section may be better understood by considering 
the zero-range limit of the confinement potential W{r) of the dot. Specifi- 
cally, let 

Vr(r) =/io(5(f) (1- (In r)fV^-), (4.35) 

in accordance with Eq. ( p.6| ). Then the Green function of H has the form 
(see (O)): 



G{r, f E) = Go(f, f E) - [Q^{E) + a^]~^G^{f, 0; ^)Go(0, r E) , (4.36) 

where = /ig ^, Go{r, r E) is given by the rhs of (|2.11|) and Qo{E) is equal 
to the expression at rhs of Eq. (|2.13|) . Hence 



Q{s- E) = Qo{E) - [Qo{E) + ao]~'Gt{s, 0; E) 



(4.37) 



and 



dQ{s- E) 
dE 



-1/2 



G,{r,s-E) 



-{QoiE) + ao)-'Go(f, 0; ^)Go(0, s; E) 
Since Go(s, 0; E) is independent of 6, we have 

'dQ{s-E)^-^'^ 



(4.38) 



Using now the results of Section 4.1, we obtain 



dQjs; E) 
dE 



dQoiE) \ , 2^ 



(go(^) + «o)"'G'o(0, s; E){Go{f, 0; E) | | Go(r-', s; E)) 



(4.39) 
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It is clear that 

(Golr-, 0; E)\Ve\ Go(f, s; E)) = Ve {Go{f, 0; E) \ Go{f, s; E)) . (4.40) 
On the other hand the scalar product {Go{r, 0; E) \ Go{r, s; E)) has the form 

{Go{f, 0; E) I Go{f, s; E)) = ! exp(-7rzeor A s) /(|f |) g{\r - s\) df , (4.41) 



and therefore it is invariant with respect to rotations of the vector s around 
the origin. Indeed, let T be such a rotation, then 

exp(— Tri^o^A Ts) /(|r |) g{\r — Ts\) dr 
= j ey.]i{—T(i^QTf /\Ts) f{\Tf\) g{\Tr — Ts\) dr 
exp(-7rz^o^A s) /(|r |) g{\r- s\) dr . 



R2 



As a result, Eqs. ( |4.39| ) and ( [4.40| ) lead to the following expression for the 
non-zero component of the Berry potential: 

dQ{s-E)\-^ fOQoiE) 



y^iP)-[^^) ^^J(<oP). (4.42) 

Using the asymptotics 

Qo{E) = 0{\n\E\), Go{s,0;E) = 0{\E\-^) as E -oo , 



we see from Eqs. ( [4.37| ) and ( ^.42| ) that in a deep zero-range well, in the sense 
that Eq <^ —JtlUc, we get 

Ve = n^0P + O{\E\^') (4.43) 
in accordance with the result ( |4.34| ) of the previous example. 



18 



5 Wilczek— Zee phase 



It was essential in the above considerations that the energy level in question 
was nondegenerate. In the opposite case the behavior of the system with 
respect to a moving perturbation is more complex, the degenerate levels 
may form different linear combinations and the change includes more than a 
simple phase factor. Nevertheless, the effect is usually labeled as the Wilczek- 
Zee phase |[WZ] . 



In magnetic systems with a homogeneous field a prime example of a 
degenerate eigenvalue are the Landau levels which constitute the spectrum 
of the unperturbed operator ( ^.11 ) with Bi = and 1^ = 0; they are 

ei= (^i + ^^huj,, i = 0,1,2,.... (5.1) 

In this section we will briefly discuss how the corresponding eigenf unctions 
behave under the influence of a moving point interaction. 

Let us first observe that the perturbation preserves the Landau levels as 
infinitely degenerate eigenvalues. Let Li be the eigenspace of H referring 
to an eigenvalue eg. It is straightforward to see that the eigenspace of Hg 
corresponding to the same eigenvalue has the following form 

Le{s) = {ijeLr. ij{s)=0}. (5.2) 

Since Li is invariant with respect to translations of the eigenfunctions, it 
is possible to select an orthonormal basis V^f ■'(s ), '?/'2^^(s ), . . . ,il)n\s), ... in 
Li{s) which depends smoothly on the point s G M^. We suppose that s is 
adiabatically moving along a smooth closed contour, s = s{t),t G [0, 1], and 
that at the initial moment t = the systems is in a state il)m{s{0)). Then 
the state 'il){t) at an instant t is given by the formula 

^(t) = e^^*/^^5^f7ia(t)^W(.-'(t)), (5.3) 

n 

where {ul^)n{t)) is a unitary matrix generalizing the Berry phase factor e*'^'-*^ 
(see ||WZ|| ). The role of Berry potential is played by the infinite self-adjoint 
matrix 

VB?)=^(^'^{^)\^A^'i'{^y). (5.4) 
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which is related to U^^\t) = U{t) by 



U-\t)U{t)) =iV},^{s{t)). (5.5) 

The solution to the equation ( [5. 5] ) along the curve C : s = s (t) is at that 
given by the path integral (the Wilson loop) 

U{C) = Vexp (^i j V{s) ds^ , (5.6) 

where V indicates a time-ordered exponential. 

The Wilczek-Zee theory has the following differential-geometric inter- 
pretation | |VDDMS|] . Consider the trivial vector bundle = x Li 



then J'f = U {{s} X Le{s) : s G M^} is a subbundle of £e with the infinite- 
dimensional typical fiber Denote by [u(oo) the Lie algebra of the unitary 
group of i"^ (the Lie algebra of skew- Hermit ian infinite-dimensional matrices). 
Then it is convenient to regard —iVmn{s) as coefficients of the differential 
form u = Ukdx^ assuming values in [u(cxd): 

CO, = {^^i\s) I v.. I ^^\s)) s = {x\x') . (5.7) 

This form is a connection form in the bundle and the operators U{C) 
are the holonomy operators in the principal W(cxD)-bundle associated with 
J^i. According to the Ambrose-Singer theorem | [KM| |, the curvature form 
VL, VL = duj + uj /\ uj determines completely the operators U{C) (the tensor 
Fjk = ifljk is the strength of the gauge potential Vk). Notice that there is an 
explicit formula (analogous to the Stokes formula) which expresses the rhs 
of Eq. (|5.6|) in terms of the coefficients of Q |[Me|| ; nevertheless, it is difficult 
to use this formula when the components of u are not commuting (which 
is the case for the matrices (|5.7| )). However, we can gain some insight into 
the behaviour of the Wilczek-Zee phase considering infinitely small loops. In 
particular, for such a loop C encircling a point sq the holonomy operator is 
given by an ordinary exponential 

U{C)=exp{nu{so)S) , (5.8) 

where S is the area encircled by the loop S. 

In the following we shall consider for simplicity the lowest Landau level 
Eq and drop the superscript for the notations. Normalized eigenfunctions 
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of the ground state Lq may be chosen in the form [CI] 

*™(n^)=(^)'^— „>o. (5.9) 

where a = sgn^o- The integral kernel Po{f,f') of the projection operator 
onto the subspace Lq equals 



Po(f, f ') = |^o|e-"^«°^^^^"e-("^-"")'/^'^o . (5.10) 
The condition ip{s) = can be then written as 

Po(s,r)^(f)df= 0, (5.11) 



and a comparison with (|5.9| ) shows that this is equivalent to 

([s,C]^o|^) = 0, (5.12) 

where [s, (] with s* G and C ^ denotes the operator of magnetic trans- 
lation ||Zal|] which acts on / G L'^{M?) as 



[s, Qfin = C exp i-m^orA s)) /(f - s) . (5.13) 

This shows that one can choose the family of the functions 

V'„^(s) = [s,l]^„, fc=l,2,... (5.14) 

for orthonormal basis in Lo{s). Let us calculate the corresponding matrix 
elements Vmn{s). It is convenient to perform the calculation in the Cartesian 
coordinates. Let r = {x,y), s= {x',y')] then 

^m(s )(f ) = exp {-m^oixy' - x'y)) ^m{x - x' , y - y') . (5.15) 

Writing as 

= ( J^V"e-"'^»''« ( (5.16) 

we find 



2™m! / V ao 



— = V& -I ,/ — Xl/ 1 — = 2— VT/ -I ,/ — VTf -, 

OX zqq flo V 2 ay 2% clq y 2 

(5.17) 
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Now we obtain from Eq. ( p.l7| ): 



dx 



— -X ,y-y) W— - a; ,y-y) 



(5.18) 



dip, 



7{s){x,y) = -Tri^oxilJmis}{x,y) + exp {-ni^oixy' - x'y)) X 



-x',y- y') W — - x' , y - y') 

zoq ao V 2 



Hence 



1 jm 



6„ 



and 



Zai aQ \ Z 



1 /T I IT \ /^e 

aQ \ Z 



2al 



(5.19) 



(5.20) 



(5.21) 



To find tlie matrix elements and ||/|^m) we make use of the 

following observations: the matrices {i{ipn \^ x' IV'm)) and («('?/'„ | Vy' I'?/'™)) are 
Hermitean, and at the same time, the numbers jxl^'m) and (\E'„ |?/|\E'm) 
are real. Taking these facts into account, we get from Eqs. ( p.20| ) and ( |5.21| ) 



% 1 Tfi 1 1 Ti 



-^{^n y "^m) + —X -Sn,m~l 

zan 0-0 V 2 



Thus 



{'^n\y\'^m) + —\I-Sm,n-1- (5.23) 



ZQn Oo V 2 



(5.24) 
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■ {y/mSn,m-l - Vn5ra,n-l) ■ (5.25) 

Since |^o|~^ = 27rao we have finally 

(^n \^x' l^m) = T^iioV'^mn + (v^f^m.n-l - Vm^n.m-l) , (5.26) 

V2ao 

(V'n |Vy l^/'m) = -7rii^oa;'5mn (v^^m.n.-l + i/m'^n.m-l) • (5.27) 

v2ao 

Because the matrices ((?/'„ |Va;' I'?/'™)) and {{ipn iVy I'?/'™)) do not commute, 
it is not easy to calculate the path integral ( p.6|) , and we turn to Eq. ( ^.8|) 
to gain some insight into the behaviour of the Wilczek-Zee phase. For this 
purpose let us calculate the curvature form VL. Since VLjk is skew-symmetric 
(w.r.t. the indices jk) it is enough to find the component VLi2. It is clear from 
Eqs. ( |5.26| ) and ( |5.27| ) that du = — vri^o'^mndx"'^ Adx^. Since the first terms in 



Eqs. (|5.26|) and (|5.27|) are scalar matrices, in order to find uj Au, we must to 



calculate the commutator of the matrices (v^ao) ^ {^/nSm,n-l — ^A^'^n,m-l) 
and —ai{^/2ao)~^ {\/n5m,n-i +v^5n,m-i) only. As a result we obtain that 

n = 27ii^Q6imSindx^ A dx^ . (5.28) 

Therefore, for an infinitely small loop C in the plane the operator U (C) 
has the diagonal matrix 

U^n{C) = diag (exp(-27rzeo5), 1, 1 , . . . , 1 , . . . ) , (5.29) 

where S is the area encircled by the loop C. Hence during an adiabatic 
evolution along the loop C the state ifjf'^ with the angular momentum m = 1 
is modified by the Berry-like factor exp(— 27ri^o'S'); the states with the others 
angular momenta m = 2, 3, . . . remain unchanged. This behaviour of the 
Wilczek-Zee phase is similar to the spectral behaviour of the Aharonov- 
Bohm Hamiltonian with an infinitely thin solenoid, which is described in 
analogy with a delta-perturbed Hamiltonian by a self-adjoint extension of a 
symmetric operator. Namely, the infinitely thin Aharonov-Bohm solenoid 
perturbs only two states with neighbour angular momenta (see, e.g., ||BV|| ). 
Similarly, in the case of the Wilczek-Zee phase the point potential changes 
two states with neighbour angular momenta: m = and m = 1. The 
opposite signs in (|0| ) and (|5.29|) can be interpreted as a "topological charge 
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conservation". More specifically, the mappings 3 i-^ [s, m = 

0. 1. . . . form a basis section of the vector bundle £q. Formula ( ^.4|) with 
m, > defines a connection in this bundle, and it is easy to show that this 
connection is flat {i.e., its curvature vanishes). Thus in accordance with the 
Ambrose-Singer theorem, all the Wilson loops (|5.6| ) are identity operators, 

1. e., the "Berry phase" for this connection is equal to zero. Adding the point 
potentials of the same strength a to each point s G we split the bundle £q 
into a sum of the line bundle Cq of the eigenfunctions in the zero-range well 
and the bundle jFg of the eigenfunctions remaining on the zeroth Landau level: 
£o = CoQ)J^o. Equations ( [4.5| ) and ( |5.29| ) show that the sum of Berry phases 
related to the summands is still zero. This effect is similar to the Berry phase 
conservation in the Born-Openheimer problem |Za2]. On the other hand, we 
have here an analogy with the Novikov formula for the Chern numbers of a 
sum of vector bundles of magneto-Bloch functions ov|| . In physical terms 
the Novikov formula states that the quantized Hall conductivity of a Bloch- 
Landau band is the sum of conductivities of the all magnetic subbands of this 
band. It remains to note that the mentioned Chern numbers are integrals of 
the curvature form. 
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